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We analyze the consequences caused by an anomalous single-particle dispersion relation suggested 
in several quantum-gravity models, upon the thermodynamics of a Bose-Einstein gas trapped in a 
generic 3-dimensional power-law potential, within the semiclassical approximation. We show that 
the critical temperature must be corrected as a consequence of such deformation in the disper- 
sion relation. Additionally, we prove that this shift in the critical temperature, in three different 
scenarios, namely, the thermodynamic limit, systems with a finite number of particles, and in a 
weakly interacting systems, can be described as a non-trivial function of the number of particles 
and the trap parameters. In other words, an appropriate manipulation of the number of particles, 
\ together with an adequate choice of the trap parameters, could be used, in principle, to enhance 

■ some quantum gravity manifestations. 

^ ■ 

I \- I. INTRODUCTION 

o 

' Several quantum gravity models suggest that the dispersion relation between energy and momentum for microscopic 
particles, must be modified as a consequence of the quantum structure of space-time 

A modified dispersion relation emerges as an adequate tool in the search for phenomenological consequences caused 
, by this type of quantum gravity models. Nevertheless, the most difhcult aspect in the search of quantum-gravity 
manifestations, is the smallness in the predicted effects Q. If this kind of deformations are characterized for some 
Planck scale, then the quantum gravity effects becomes very small [H, Q. In the non-relativistic limit, the deformed 
dispersion relation can be expressed as follows 0, Q (in units where the speed of light c — 1) 
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O ; Where Mp (~ 1.2 x lO^^eV) is the Planck mass, m is the mass particle in question, and p is the momentum. The 
three parameters ^i, ^2, and ^3, are model dependent [Hill, and should take positive or negative values close to 1. 
Equation ([T|) is the starting point in the search for small manifestations of quantum gravity at low energies. 
^ ' The use of Bose-Einstein condensation phenomenon as a tool in the search for some quantum-gravity manifestations 
in our low energy world, is not a new topic, and has produced an interesting series of studies in this direction i3- 
However, as mentioned above, the corrections caused by some quantum-gravity manifestation, in the relevant 
properties associated to the Bose-Einstein condensate, are still very small to be measured, then the main goal is 
search the possibility to enhance such possible manifestations. 

The analysis of a Bose-Einstein condensation phenomena in the ideal case, weakly interacting, and with a finite 
number of particles, trapped in different types of potentials [T2| - [23| shows that the main properties associated to 
the condensate, and in particular the critical temperature, depends strongly of the characteristics of the trapping 
potential in question, together with a non-trivial function of the number of particles. A logical step, is introduce a 
modified bosonic gas in this scenario and analyze the behavior of such systems in this quite general description. For 
this propose, let us define the next "modified Hamiltonian" 



H = ^ + ap + U{r). (2) 
2 m 

Where p is the momentum, m is the mass of the particle and the term ap is the leading order modification in 
expression ([I]). 

The potential term 
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in the Hamiltonian ^ is the so-called generic 3-dimensional power-law potential, where Ai and are energy and 
length scales associated to the trap [23j . On the other hand, are the d radial coordinates in the rii-dimensional 
subspace of the 3-dimensional space. The sub-dimensions rii satisfy the following expression in three spatial dimen- 
sions 



d 

J2n,^3. (4) 

If in equation ^ d = 3, ni = 122 = = 1 then the potential becomes in the so called cartesian trap, if d = 2, 
ni — 2 and ^2 = 1 then we obtain the cylindrical trap, and if d = 1, ni = 3, we have the spherical trap. On the 
other hand, if s; — >■ 00 we have a free gas in a box. In this sense the potential included in the Hamiltonian ([2|) is quite 
general. Different combinations of these parameters gives different classes of potentials according to ([3]). 

The outline of the present paper is the following one. In section 2 we deduce the shift in the critical temperature in 
the thermodynamic limit, in function of the number of particles. Section 3 studies the shift in the critical temperature 
caused by a finite number of particles. Section 4 analyzes the properties of a weakly interacting Bose-Einstein gas, 
for systems with a large number of particles. Finally, section 5 analyzes the main results and adds some comments 
concerning the perspectives of precision tests with Bose-Einstein condensates. 



II. DEFORMED BOSONIC GAS TRAPPED IN A GENERIC 3-DIM POWER LAW POTENTIAL IN 

THE THERMODYNAMIC LIMIT 

Let us calculate the number of microstates for the modified Hamiltonian In the semiclassical approximation 
the number of microstates [25j , and the semiclassical energy associated with the modified Hamiltonian ^ are given 
respectively by 



+ ap + U{r). (6) 

Zm 

From equation ([6]) we obtain the relation between energy and momentum, 



p = (2m)i/2 + Irna^ _ ^ ^ am. (7) 

i—l ' 



Integrating expression ((S]) in the momentum space, using the volume of a sphere, together with equation (O at first 
order in a leads to 



4-Tr , /■ / 1 \ 3/2 

— aj^-^^m{2m) J d'^ r(^e + —ma^ — U {r)J . 

The integration over the space in equation ^ using the definition of the Beta function B{a,b) — a;°^^(l — 
x)^~^dx = r(a)r(6)/r(a + b), after some algebraic manipulation, leads to the number of microstates associated to 
this system 
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Where C is a constant that depends of the potential in question (in the case of a harmonic oscihator C = 8), and 
r(a;) is the Gamma function. 

The number of microstates per energy unit can be obtained by taking the derivative with respect to e in equation 

m 



where 



n{e;a) = — ^-^=fe+-ma^j Qi - aU + -ma'' ] n^. 
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Let us calculate the transition temperature. The number of particles in the continuum approximation is given by 
the following expression in the grand canonical ensemble (2^ 



(13) 



where No is the number of particles in the ground state, and z is the fugacity which is related to the chemical 
potential through the expression z = exp(/3/i), here /3 = 1/kT, and k is the Boltzmann constant. 

Using equation ((TU)) and equation p^ . assuming that mo? 12 << kT, with the change of variables x = (3e = 

/3(^e + ^ma^^ , allows us to express the number of particles as follows 



N = No + CntiA; aFT^^ + 1 
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(15) 



is the parameter that defines the shape of the potential, being Ze// = z exp /3{ma'^ /2) an effective fugacity. Addi- 
tionally T{x) is the gamma function. 

The functions g^{z) are the so-called Bose-Einstein functions defined by [2^. l25j 



r(i/) Jo z-^e^ - 1 



(16) 



The Bose-Einstein function (|16p diverges for z = 1 when i/ < 1 [241] . 

Expanding around a = to first in order a, using the properties of the Bose-Einstein functions ^4 
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allows us to write expression ([14]) as 



si 



(18) 



Let us define 



si 
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as the characteristic volume associated with the system . We can see that from expression ([T^ . when -> 00, then 
Vchar becomes the volume associated with a free gas in a box. In this sense Vchar can be interpreted as the available 
volume occupied by the gas [3, HI] . 

When the number of particles, N ^ is large, we may neglect the zero-point energy contributions and we can equate 
the lowest energy to zero [26j . Hence, for the onset of the condensation in the thermodynamic limit the number of 
particles in the ground state is negligible TVq = and the fugacity z = 1, with these assumptions the Bose-Einstein 
functions g-yiz) are given by the Riemann Zeta functions ^(7) [131 • is noteworthy to mention that the phenomenon 
of condensation for a modified bosonic gas in the thermodynamic limit is possible when 7 > 3/2 according to the 
properties of the Bose-Einstein functions (|16p. Expression (jl8l) at the transition temperature now becomes in 
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where T^, is the critical temperature. 

If we set a = in equation (pn)) we obtain the usual result [2E 



N = Vchar 



(21) 



From which we obtain the usual definition of the critical temperature Tq in the thermodynamic limit 
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At this point it is noteworthy to mention that the most general definition of thermodynamic limit can be expressed 



iV ^ 00, Vchar ^ 0, (23) 

keeping the product NVchar — ^ const, and is valid for all power law potentials in any spatial dimensionality f2lj. 
Using these facts, the critical temperature in the thermodynamic limit is well defined. 
Equations (1^ and ([^ implies 
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It is noteworthy to mention that the correction in the critical temperature (j24p depends strongly on the functional 
form between the number of particles and the parameters of the potential. 
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Let us suppose now that our condensate is trapped in an anisotropic three-dimensional harmonic-oscillator potential. 
For this trap, the value of the shape parameter is given by 7 = 3, with — hu!i/2, = y^h/mLUi, (see expression 
(fT9| ). and with the definition ui = (a;ia;2W3)^/^, we obtain from expression ([24| 
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If we set a = in equation then there is not shift in the critical temperature, = 0. In typical experiments 
the number of particles vary from a few thousand to several millions, and frequencies Cojlix from tens to hundreds 
of Hertz Q. For instance, in the case of with say, wi = 27r x 50Hz, W2 = 27r x 40IIz, ^3 = 27r x 20IIz, and 

a = ^i-^-c > 0, c is the speed of light and Mp is the Planck's mass, with say, N = 10^^ and N — 10^ we obtain from 

([25| that the shift in the critical temperature is given respectively by ~ 'fi9.7 x 10^^ and ~ i^i3 x 10~^. 



III. DEFORMED BOSONIC GAS AND FINITE NUMBER OF PARTICLES 

In this section let us calculate the leading correction in the critical temperature caused for a finite number of 
particles in our modified bosonic gas. For this purpose we must take into account the associated zero-point energy 
eo [231, . At the critical temperature the chemical potential is given by /x = eo |26|] . Let us start with equation p4)) 



TV = A^o + Vchar 
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where z^ff = e'^(^+™"^/^' is an effective fugacity, and gy{x) is the Bose-Einstein function. 
We can expand the last equation using the behavior of Bose-Einstein functions about zero [24 
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In expressions (l27l) and (|28|) . r(a;) is the gamma function and C,(y) is the Riemann Zeta function. Expresion ((27|) is 
valid for all < 1, and to all non-integral v > \, and expression (1^51) is valid for positive integers m, both expressions 
are valid in the limit 5 — s> 0. The case = 1 is rather simple, in this case we have 



5i(e-^)--ln(l-e-*)^ln(l/(5), 



(29) 



when (5-^0. 

With this assumptions, we can express equation (|26|) at Tc, with Ze//e = 6^"'''=) as follows (at Tc the chemical 
potential is bounded i.e., /i = eo) 



/ m 
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Equation (|30)) implies 
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Using the result for the critical temperature Tq in the thermodynamic limit expression (1221) and equation (j3ip . we 
obtain 
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(32) 



If we set a = in equation ([5^ we recover the "long range effects" corrections given in [23| . For an anisotropic 
three-dimensional harmonic oscillator potential (7 = 8 and using the expression ([32]) together with expressions ((27)) 
and (|28p . we obtain that the relative correction in the critical temperature first order in a and eo, in function of the 
number of particles is given approximately by 
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If in equation ([55)) we set a = 0, we obtain the first correction caused by a finite number of particles [I3lll5lll6ll26l|. 

In fact, the term proportional to the zero-point energy eo, is the leading correction caused by the finite size effects of 
the system. 

For instance, in the case of the ®''i?5, with a — Ci^IF"^ > 0, c is the speed of light and Mp is the Planck's mass, 
we obtain that the leading correction in the shift ([55)) caused by the deformation for N = 10^ particles, is given by 
(4.5 X 10-7(1)1/2) + 3.7 X 10-20(-^)) , and for N = 10^ we obtain ^ (l-4 x 10-6(1)1/2) + i.n x io-i6( Jfl^)^ , 

where luq = |(wi + W2 + ^^3)- 

This correction, unfortunately is small and show the requirement not only of a finite system, but also with a small 
number of particles, in which the approach followed in this paper is an approximation, due to problems of fluctuations 
and non-equivalence of the ensembles when the system contain both, a finite and small number of particles 0, Q . 



IV. WEAKLY INTERACTING MODIFIED BOSONIC GAS 



A realistic system does include interaction among its constituents. In this aim let us estimate the correction in 
the critical temperature for a deformed bosonic gas taking into account interactions among its constituents within 
the Hartree-Fock theory in the semi-classical approximation, which basically consist in the assumption that the 
constituents of the gas behave like a non-interacting bosons, moving in a self-consistent mean field and is valid when 
e{r,p) » /i for dilute gases [sl. [26|. 

Let us propose a particularly simple modified Hartree-Fock type spectrum. 



2 

e(f,p) = ^ + ap + U{r) + 2U^n{r). (34) 

Where the term 2Uon{r) is a mean field generated by the interactions with the other constituents of the bosonic 
gas Q. C/q is the coupling constant which is related with the s-wave scattering length a through the next expression 

47rft2 

Uo = a. 35 

m 

On the other hand n{r) is the spatial density of the cloud. 

Let us calculate the spatial density associated with our modified Hartree-Fock spectrum (15^ . In the semiclassical 
approximation, the single-particle phase-space distribution is given by [1, [2^ 

The number of particles in the 3-dimensional space obey the normalization condition. 



7 



N = - — — - / d?rd?pn{r,p) = / (frn{r) = / d^pn{p), (37) 



where 



n{r) = I (fipn{f,p), n{p) ~ / (firn{r,p) 



(38) 

are the spatial and momentum densities, respectively. 

Expression p7p is equivalent to expression ()13p . In expression (j37]) we assume that the number of particles in the 
ground state above the critical temperature is negligible. Using the expression (|34|) . and integrating expression (|36|) 
over the momentum space, allows us to obtain the spatial distribution associated with the modified Hartree-Fock 
spectrum (p4|) 
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is the de Broglie thermal wavelength, /Xejj = A* + and 5t/(a;) is the Bose-Einstein function. If we set a = 

in equation (1391) we recover the usual result for the spatial denstity, in the semiclassical approximation [s], [2^ . 

Using the properties of the Bose-Einstein functions, equation (jl7p . we can expand expression (j39p around — 0, 
with the result 



n{r) « no(f) + Uo{2kT)-^\-'' [g3/2(^)5i/2(^)] + aC/o(2KT)-iA-5 ^g,/^{Z)go{Z) + gi{Z)gy2{Z)\ . (41) 

where 



Being no(r) the space density distribution for the ideal case Uq = 0, 



(42) 



(43) 



Integrating the normalization condition (1371) using expression (|4ip with the corresponding potential allows us 
to obtain an expression for the number of particles in function of the chemical potential ji, the temperature T, the 
coupling constant J7o, and the deformation parameter a, given by 
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Where 
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Where, Ze// = e^^'^^"^° If we set ?7o = in equation (|44p . we recover our previous result for a modified ideal 
bosonic gas in the thermodynamic limit (jl4p . Additionally, setting a = and Uq — 0, we recover the usual result for 
for an ideal bosonic gas in the thermodynamic limit ([21]) . 

In order to obtain the leading correction in the shift for the critical temperature caused by the interactions in our 
deformed bosonic gas, let us expand the expression (j44p at first order in T = Tq, ~ 0, Uq = 0, and a = 0, Tq is the 
critical temperature in the thermodynamic limit given by expression p2|) . with the result 
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At the critical temperature for large A^, in the mean field approach fic — 2 UonQ^r — 0) which means that the 
critical density at the center of the trap is the same as for the non-interacting model. In the usual case when a = 0, 
no{r ^ 0) = A^'^C(3/2) in the large A^ limit 5], but in our case we have to modified the value for /ic at the critical 
temperature according to expression (|43l) , because of the divergent behavior of the Bose-Einstein functions related to 
no(r = 0), namely gi{Z) and 51/2(2') when Z = 1. When the integrals associated with the Bose-Einstein functions 
converges, the value mo? jl is negligible and can be replaced by zero. However when the integral associated to the 
Bose-Einstein functions can diverge at Z ^ 1 the minimal energy associated with the corresponding system must be 
taken into account [21j. 

In order to avoid the divergent behavior of these functions, we have to take into account the minimal energy of this 
system. We are interested in the corrections due to a in the large A^ limit, so we will take as a minimum energy of 
the system ma^/2, ignoring the corrections due to the corresponding ground state energy or finite-size effects, which 
have already been calculated in equation ([32]) . 

Let us define no(r = 0) at the critical temperature using expression (j43p as 



no(f = 0) = A-353/2(e^^'""'/^) - 2af/oA-^ (^).gi(e'^=-^/2). 



(47) 



With the help of expressions (P7|) . and (23 |. we can define the Bose-Einstein functions 53/2(6'^^™" and 
5i(e'^'='""'/2), ^j^gj^ (/3^maV2) ^ as 



ff3/2(e^=™"^/^) ^ C(3/2)+r(-l/2)(^) 
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(48) 
(49) 



Neglecting the second order terms in Uq and a, allows us to write fic using expression (|47l) as follows 



(50) 



If in equation (jSOp we take the limit when a ^ we recover the usual value for /i^ at the critical temperature @ . 
Inserting ()50p in (|46p . we finally obtain that the shift in the critical temperature caused by the interactions in a 
modified bosonic gas in function of the number of particles is given approximately by 
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Setting a = in equation ([ST|) we obtain the usual shift in the critical temperature caused by weakly interactions. 

Let us analyze the expression ([STj) . The first right-hand side term is associated to the usual correction caused by 
interactions, the second term is a correction caused by the deformation term, and the last one, corresponds a correction 
caused by interactions and the deformation term. A decreasing or increasing in the temperature depends of the sign 
of a and a. For a > (repulsive interactions) and a > 0, the correction in the shift caused by the deformation is a 
positive shift. On the other hand for a < (attractive interactions) and a > the shift caused by the deformation 
term is negative. In the case of a > and a < we obtain a negative shift, and for a < and a < 0, the second 
term is negative and the third is positive. An increasing or decreasing corrections, depends of the sign of a and a 
and its combinations. The second term in expression (j51|) proportional to , becomes smaller when the number 

of particles is large but, the third term proportional to the scattering length becomes bigger. 

In the case of an anisotropic three-dimensional harmonic oscillator potential 7 = 3, with Ai = huJi/2, ai — ^ h/muii, 
then from expression (jSip we obtain 
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Where we used the usual definitions 
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By analyzing expression ([5^ . we obtain for with say N — 10^*, N — 10^, and 10^ particles, the correc- 

tions caused by the deformation term are given respectively by Ci(^3T724.5 x 10^^ + 9.12 x 10^^1n(3.7 x 10^^ l?)^, 

(3172 1-42 X 10-^ + 9.12 x 10-9ln(3.7 x 10"^)), ^ (-1/24.5 x 10"^ + 9.12 x 10-^ln(3.7 x 10i"|-)), with the 

scattering length a — 5.77 x 10^^ m and ui = (a;iaj2W3)^/'^. 

On the other hand, with say, ui Ki 2tt x 300 Hz, for the third term proportional to the product aa becomes 

bigger than the first one proportional to a in equation (I53p . when the number of particles N is approximately iV > 10^, 
which means a change of sign when a < 0, assuming ^1 ?» —1. 

If we set a = in expression (|53p we recover the result given in [l^ . 



V. CONCLUSIONS 



Using the formalism of the semiclassical approximation, we have analyzed the Bose-Einstein condensation for a 
modified bosonic gas trapped in a 3-D power law potential in three regimes, namely in the thermodynamic limit, in 
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systems with a finite number of particles, and in a weakly interacting systems in the large N limit. We deduced the 
density of states for the aforementioned system in the ideal case and we have deduced the parameter that characterizes 
the shape of the potential, which are directly related to the existence of the condensate in the thermodynamic limit. 
We also deduced the shift in the critical temperature in the thermodynamic limit, in finite systems, and for a weakly 
interacting bosonic gases in function of the number of particles, equations (|24l) . (|32l) . and ([STj). which are valid for 
any potential defined by the generic 3-dimensional power-law potential ([3]) within the semiclassical approximation. 
Nevertheless, the case of an ideal Bose-Einstein gas in a box in three dimensions, deserves a separately analysis. 

In the case of an uniform Bose-Einstein gas in a box Si — > oo, and hence 7 = 3/2, then, according to equation (|18p 
the number of particles for this system is given by 



N = No+ (^KTy^^fligy^iz) - a(KT)n2gi{z) + a"" (kT^ ^^g^/^iz) (55) 

Apparently Bose-Einstein condensation is not possible for this kind of systems in three dimensions in the thermo- 
dynamic limit, because of the divergent behavior of Bose-Einstein functions 51(2:) and gi/2{z), when /i = at the 
transition temperature [ll| . In the usual case, when a — Q, we have only the first term in expression (|55l) . in this 
situation the condensation is possible in three dimensions for a bosonic gas in a box, and assuming that the minimal 
energy for this system is zero, we can fix the value of /i = safely without problems of divergence and (73/2(1) is just 
C(3/2). Nevertheless, the minimal value of the energy for a modified bosonic gas in a box is not zero, in fact, this 
value is mo? 12. 

These assumptions allows us to write expression ([55)1 at the critical temperature as follows 

N ^(KT})'\^g^l2{ef''^'^"^'')-c.(KT^^2gx{^^^^ (56) 
With the help of expressions we may now express equation (15^ as follows 



N = (AtT,)3/2iiiC(3/2) + a{KT,) 



r]ir(-i/2)mi/^ 
21/2 



^2 ln(- 



2Krc\ 2i/2f23r(i/2) 



V ma 



7,1/2 



(57) 



If we set a = in expression (|57|) . we obtain the usual result for a bosonic gas in a box 24]. Using the definition of 

the critical temperature for a bosonic gas in a box in three dimensions Tq — 2vmK ( ^^3/2 ) i-"^ thermodynamic 

limit [13, [IB] and with the expression (|57|) . we obtain the shift in the critical temperature caused by the deformation 
term in function of the number of particles 
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Where V is the volume of the container. Expression ([58)) shows that the condensation for a modified bosonic gas 
in a box is possible, taking into account the minimal energy of the system in question. Nevertheless, this system is 
unstable, because the fluctuations in the number of particles becomes anomalous [T]| . 

In the case of the different values for the shape parameter, we notice that for spherical traps d = 1, ni = 3, and 
setting, say si = 2 then, 7 = 3, hence from equation (|24p in the thermodynamic limit we obtain 



^ - aN-'/\ (59) 

On the other hand for cylindrical traps d = 2, ni = 2, and 71.2 — 1, setting, say Si = 2, and S2 = 1, the trap 
parameters are given respectively by 7 = 9/2. From equation (j24p in the thermodynamic limit we obtain 



(60) 
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These two simple cases, illustrate how the tiny contributions caused by the quantum-structure of the space-time, 
could be enhanced, varying the number of particles. 

A very important consequence, is that for values between 3/2 < 7 < 5/2 we have condensation, but the system 
is unstable |11| . In this range of values of the shape parameter 7 we have, for example, potentials of the type 
V{f) ^ + + in the case of a Cartesian potential, V{f) ^ + z'^ in the case of cylindrical traps, and V^(r) ~ 
in the case of spherical traps. We have noticed also, that for potentials of the type ^(r) ^ p"' + z^, V{r) ^ or any 
other combination such that 3/2<7<5/2we have condensation, but the system is apparently unstable (assuming 
that the zero point energy is negligible), in the thermodynamic limit, because the fluctuations in the number of 
particles becomes anomalous [11|. When the fluctuations in the number of particles becomes anomalous, in the 
thermodynamic limit, could leads to an immediate collapse or implosion of the condensate [2lj. This last assertion, 
allows us, in principle, obtain a criterium to discriminate when the quantum-gravity manifestations are more feasible 
to be measured, by analyzing the stability for these systems [lH. On the other hand, we prove that, in the case 
of stable systems (for example, the harmonic oscillator) , the tiny contribution due to the quantum-structure of the 
space-time is enhanced by the number of particles N. 

The present work shows that the effects caused for the quantum structure of space-time could be, in principle, am- 
plified taking into account corrections caused by a finite number of particles in an ideal bosonic gas, or a large number 
of particles in an interacting bosonic gas, together with an adequate election of the parameters in the corresponding 
potential. The shift in the critical temperature of a condensate could provide bounds for the parameters associated 
to the deformation suggested in several quantum gravity models. 
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